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Abstract 

We present first-principles studies on the orbital states of the layered perovskites Ca2_xSr2,.Ru04. 
The crossover from antiferromagnetic (AF) Mott insulator for x < 0.2 to nearly ferromagnetic (FM) 
metal at x = 0.5 is characterized by the systematic change of the xy orbital occupation. For the 
AF side (x < 0.2), we present firm evidence for the xy ferro-orbital ordering. It is found that the 
degeneracy of t2g (or eg) states is lifted robustly due to the two-dimensional (2D) crystal-structure, 
even without the Jahn- Teller distortion of RuOg. This effect dominates, and the cooperative 
occupation of xy orbital is concluded. In contrast to recent proposals, the resulting electronic 
structure explains well both the observed X-ray absorption spectra and the double peak structure 
of optical conductivity. For the FM side {x = 0.5), however, the xy orbital with half filling opens a 
pseudo-gap in the FM state and contributes to the spin 5'=l/2 moment (rather than 5=1 for x=0.0 
case) dominantly, while yz, zx states are itinerant with very small spin polarization, explaining the 
recent neutron data consistently. 

PACS numbers: 71.27.+a, 71.70.-d, 74.70.Pq 
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I. INTRODUCTION 



The spin and orbital structures in Ca2-xSra;Ru04 (0< x <0.5) attract much attention 
recently as a model system of electrons. It is an antiferromagnetic (AF) Mott insulator 
for X <0.2, but nearly ferromagnetic (FM) metal at x=0.5^^^. There are four electrons 
occupying t2g states per Ru^+, while eg states are empty. The isovalent substitution of Sr 
for Ca does not change the number of electrons, but modifies the crystal structure systemat- 
ically^- due to their different ionic radii. The key issue here is the orbital degree of freedom, 
which couples with lattice and magnetism strongly. Extensive experimental^'^'^^ and theo- 
retical^iiSiiLi^ studies have been done for this issue, while the present understanding of this 
compound is quite controversial. 

First, for the AF side (a;=0.0), the preferential occupation of the xy orbital was pre- 
dicted by earlier calculations^iiS. However the O Is X-ray absorption spectroscopy (XAS) 
study reported the occupation of about 0.5 holes in the xy orbital^. In addition to this, 
the photoemission experiment suggested strong spin-orbit coupling (SOC), leading to the 
proposal of complex orbitala^. The presence of 0.5 holes in the xy orbital also led Hotta 
et al.^sM to propose the "antiferro-orbital ordering" (AFO) with 2x2 periodicity of RuOg in 
the plane. The recent experiment on optical conductivity was interpreted in terms of this 
orbital pattern^. However, the lattice distortion corresponding to the AFO (i.e., the 2x2 
structure) is not observed experimentally, and the orbital pattern of AFO is not consistent 
with the recent resonant X-ray scattering (RXS) experiment)^, in which no superlattice peak 
was observed. Second, for the nearly FM side (x=0.5), the fitting of the Curie- Weiss type 
susceptibility suggested 5=1/2 spin moment^, in sharp contrast with the high-spin (5*=!) 
state for x=0.0. A similar picture may hold also for a;=0.2''. The recent polarized neutron 
diffraction studies for x=0.5^ suggested dominant magnetization distribution on the xy or- 
bital rather than the theoretically suggested yz/zx orbitals^*^. Therefore it is an important 
and challenging issue to determine the orbital state of Ca2-xSr2;Ru04 explaining consistently 
the available experiments, which we undertake in this paper. 

We show that the crossover of magnetic properties with increasing doping x is charac- 
terized by the systematic change in the xy orbital occupation. For the AF side (x <0.2), 
we conclude that the predicted ferro-orbital (FO) ordering with dominant xy occupation 
explains well (a) the XAS spectra, and (b) the double peak structure of the optical conduc- 



2 



tivity. We further point out that the stabihzation of the xy state is quite robust, and mostly 
due to the two-dimensional (2D) crystal-field, which has not been considered seriously so 
far. In contrast to the AF side, for x=0.5, the electronic structure is characterized with the 
yzj zx states with small moment, and the xy orbital with half filling, which opens a pseudo- 
gap in the FM state and contributes to the S'=l/2 moment dominantly. This picture, which 
is in strong contradiction to the previous proposaii^, explains the neutron data^ consistently. 

II. CALCULATION METHOD 

The calculations were done with the first-principles plane-wave pseudopotential method. 
In our previous work™, the effects of structural distortions were emphasized by taking into 
account important distortions mostly, however here we use the experimental structures to 
compare the calculated results with the measured spectra. Five experimental structures are 
considered, i.e., x=0.0 at IIK and 295K, x=0.1, 0.2 and 0.5 at lOK^. Two kinds of mag- 
netic states, FM and staggered AF states, are considered for each structure. In the previous 
work^, the local spin density approximation (LSDA) was adopted for the entire range. The 
validity of LSDA for predicting the stable magnetic state of ruthenates was demonstratecii^. 
However for Ca2Ru04, though LSDA can predict the stability of AF magnetic order, the 
band gap does not oper^*li. As we aim to analyze the optical conductivity also, a proper re- 
production of the band gap is important. Therefore, we employ the LDA-|-[/ schemei^>i^ for 
the insulating region (0< x <0.2) with [/eff=2.5eV to reproduce the band gap of Ca2Ru04. 
The LSDA was employed for other region of x. The details of other aspects of the calcu- 
lation were described in the previous paper—. For the studies of XAS spectra and optical 
conductivity, the interband transitions are calculated from the converged Kohn-Sham wave 
functions and eigen values, by using the core compensation form of pseudo-wave functions 
(see Ref.— for technical details). 

III. RESULTS AND DISCUSSION 
A. Ca2Ru04 (x = 0) 

We start with Ca2Ru04 at lOK. Our calculations predict the AF ground state with 
energy gain about 28.2meV with respect to the FM state. The mechanism for stabilizing 
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the AF state was analyzed in our previous paper— in terms of lattice distortion, flattening, 
rotation and tilting of RuOe- The projected density of states (PDOS) in Fig. 1(a) suggests 
that the present situation may correspond to the localized spin picture of 5*=!, where the 
xy states are fully occupied and the yz/zx orbitals are half filled. In reality, as the Ru 
Ad orbitals are extended, inter-site hybridization modifies the picture quantitatively. For 
example, the occupation numbers (in LDA+U) are given as follows: nj.j^=0.86, n'yz/zx~^-^'^ ^ 
n^j^=0.79 and nyz/zx~^-'^^^ with up and down arrows indicating the spins (see Figs. 1(a) and 
2(a)). The cooperative occupation of the xy orbital for all Ru sites forms the FO ordering. 
These occupation numbers give 1.25/1^ as the magnetic moment of a Ru atom Mpiu_totai, 
which agrees well with the experimental value of l.SfiB^- As for the spin polarization of 
oxygen, the in-plane oxygen atoms have no net polarization due to symmetry while the apical 
oxygen atoms are strongly polarized (about 0.1 muB) due to the strong p — d hybridization 
because the yz / zx orbitals which extend toward apical oxygen contribute dominantly to the 
Ru spin polarization. The corresponding values in LSDA for the above orbital population 
are n|,j^=0.80, nlyz/zx~^-^'^^ ^iy^^-'^^ ^^^^^^=0.40. Therefore, the basic feature of the 
orbital population remains the same in LSDA, though the orbital polarization and the 
spin polarization are reduced compared with those in LDA+f/. For example, LSDA gives 
0.93/iB as MRu_totai- The reduction in these polarizations is due to the stronger inter-site 
hybridization caused by the vanishing band gap in LSDA. Therefore, we believe that the 
results obtained by LDA+U which adjusts the effective Coulomb repulsion Ucs to reproduce 
the band gap should be quantitatively more reliable than those by LSDA. The present results 
for the electronic structure are consistent with previous calculation o^°i^^ . 

As for the orbital character, there are three factors mainly contributing to the stabilization 
of the xy orbital: 

1) The energy level splitting due to 2D crystal field in the layered structures (in contrast to 
its three-dimensional (3D) counterpart, where the t2g (or eg) states are degenerate). Fig. 3 
(main panel) shows the calculated PDOS for the non-magnetic (NM) state of Ca2Ru04 
without the Jahn- Teller distortion (i.e. three Ru-0 bonds having equal length). In the 
inset of Fig. 3, we show the calculated energy difference between the xy and the yz, zx band 
as a function of the tetragonal JT distortion, by calculating the center of gravity for each 
band. Clearly, even for the hypothetical tetragonal Ca2Ru04 without JT distortion, the xy 
orbital is lower in energy than the yz and zx orbitals by about 0.2eV as estimated by our 
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calculation. This lowering of xy state energy is added to that due to the tetragonal distortion 
as demonstrated in the almost uniform downward shift of the dashed line (Ca2Ru04) from 
the solid line (CaRuOa) in the inset of Fig. 3. As a first approximation, the origin of the 
energy level splitting is attributed to the following geometrical aspect in the second neighbor 
configurations: in Ca2Ru04, Ru-O-Ru in the ah plane is replaced with Ru-O-Ca along the 
c axis. Two factors coming from this geometrical aspect contribute to the energy level 
splitting. First, the electrostatic potential due to the reduced positive charge of the second 
neighbor Ca^+ along the c axis will raise the energies of those orbitals extending along the 
c axis, like the yz and zx (or — r^) states. This can be seen in Fig. 3 from the global 
upward shift of yz, zx bands and their oxygen bonding counterparts relative to those of the 
xy state. Second, due to the absence of Ru-0 bond at one side of the apical oxygen along 
the c axis, the remaining Ru-0 bonds are strengthened. This will further push up yz and 
zx states, which are the anti-bonding parts of the 2p-4:d hybridization. 

2) The compressive JT distortion. The xy state is further lowered by about 0.06eV (see 
inset of Fig. 3) due to the 2% shrinkage of the apical Ru-0 bond length observed for Ca2Ru04 
at lOK. 

3) The orbital-dependent hybridization. In the AF state, the inter-site hybridization 
between the occupied orbitals and the unoccupied ones, which is the origin of the super- 
exchange, will push up the unoccupied yz/ zx states, and again enhance the splitting between 
xy and yz/zx in minority spin by about 0.07eV from our calculations. This effect does not 
exist in the FM state, and will explain the reduced occupation of xy orbital in the FM 
solution even for x=0.0 (see Fig. 1(b)). 

As the results, nl.y reads 0.79, 0.72, and 0.67 for the cases corresponding to (x=0.0, 
IIK, AF), (x=0.1, lOK, AF: JT distortion is nearly vanishing), and (a;=0.0, IIK, FM), 
respectively. It is clear that the xy occupation is dominantly determined by the energy level 
splitting due to the 2D structure, which is comparable with the typical band width of t2g 
states. The existence of such 2D crystal-field is common for all the layered perovskites, while 
it has not been considered seriously so far. This effect, which is not taken into account in 
the analysis by Hotta and Dagottc^^, will certainly suppress the stability of their proposed 
AFO state for Ca2Ru04. 

We also performed self-consistent calculations by including the SOC using the relativistic 
fully separable pseudopotentialai^ in the framework of non-coUinear magnetism. We found 
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that the obtained orbital occupations are almost identical to those for the case without 
SOC, in the AF gound state of Ca2Ru04. This is consistent with above discussions in the 
sense that the crystal field is strong enough to quench the orbital moment and stabilize 
the real-orbital xy. Note that the photoemission measurement^ was done at 150K which is 
above the Neel temperature where the degree of orbital polarization is reduced. 

The electronic structures obtained in the present work can account well also for the 
observed XAS^ for Ca2Ru04 though the number of holes in the xy state is not exactly 0.5 
(the estimated number of holes is between 0.2 and 0.3 from our calculations for the AF 
ground state). There are two important aspects in the experimental analysis. One is the 
distribution of holes among t2g states and the other is its variation from the low-temperature 
AF phase to the high-temperature paramagnetic (PM) phase. Such information was derived 
from the observed relative intensity of two XAS peaks (A and B) as a function of light 
incidence angle 9. The peaks A and B correspond to the ls-2p transition at the apical and 
the in-plane oxygen, respectively. In the single electron transition limit, we obtain the XAS 
spectra from the converged Kohn-Sham eigen states by calculating the matrix elements of the 
optical transitioni^. For the low-temperature AF phase, the calculation is straightforward. 
However, we have to introduce an approximate treatment for simulating the XAS in the 
high-temperature PM state. Here we use the following simplest approximation in the cluster 
expansion technique for treating the alloy problen>i^, assuming that a magnetic state can 
be regarded as an alloy with two constituents, up spin state and down spin state. As the 
simplest approximation, we use a dimer for the cluster. Then the XAS of paramagnetic 
state can be approximated simply as an average of XAS of FM and AF states. (If we use 
a larger cluster, other magnetic orders have to be taken into account and the short-range 
order effect can be included.) As shown in Figs. 4(a) and (c), our calculated XAS at IIK and 
295K can be well compared with the experimental spectra taken at 90K and 300K (Fig. 3 
of Ref.— ). At elevated temperature, which tends to suppress the compressive JT distortion 
and enhance the spin disordering (increasing FM component in spin configurations), the n\.y 
will decrease as obtained in our calculations for the 295K structure (see Figs. 1(c) and (d)), 
being consistent with the experimental tendency. 

Recently, the anisotropic optical conductivities in Ca2Ru04 were measured by Lee et al^ 
and Jung et alM. For the E//a spectra (see Fig.5), two peak structures (called a and /? 
at about 1.0 eV and 2.0 eV respectively) are related to the 4d — 4d transitions (the strong 
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peak at 3.0eV is due to the p — d charge transfer). By increasing the temperature, spectral 
weight transfer from the strong peak j3 to the weak peak a is observed and simultaneously 
the gap is reduced. We calculatei^ the optical conductivities by using the Kubo formula^i, 
and all transition matrix elements are calculated from first-principles^^. As shown in Fig. 5, 
the experimental results can be well explained by our calculations, and the corresponding 
transition paths are indicated in Fig.l by red arrows. Again we use the simple average of 
FM and AF solutions to simulate the PM state at 295K. For the low temperature AF state, 
the peak (3 at 2.0 eV is mostly due to the transition from the occupied majority spin yz (zx) 
state at one Ru'^"'" site to the unoccupied minority spin yz (zx) state at the neighboring Ru^+ 
site. In the simple energy level picture (i.e, the occupation number of each orbital is regarded 
as integers, either or 1), it is easy to understand the origin of peak (3 (which corresponds to 
transition shown in Fig. 1(c) of Ref.— with energy U + J), but not that of a. This argument 
has been used by Lee et al.— against the FO ordering. So they pursue the proposal for the 
AFO stat©'''^, which has been shown in the above discussions not to be favored. However, 
our calculation clearly shows the existence of peak a at about l.OeV, which is due to the 
no n- vanishing transition matrix element between xy and yz/ zx states, and crucially depends 
on the orthorhombic distortion and the admixture of xy and yz/ zx states of Ru-4(i. One of 
the strong contributions to a is the transition from the occupied majority spin xy state at 
one Ru^+ site to the unoccupied minority spin xy component mixed into the yz^ zx states 
at the neighboring Ru''"'' site (see Fig. 1(a)). Therefore, the double peak structure is just 
the natural consequence of strong hybridization. The peak (3 should be strong at the low- 
temperature AF phase, and sensitive to the magnetic ordering. Actually, this transition 
should be suppressed in the FM configuration (as shown in Fig. 1(d)), because of the flipping 
of the spin. The peak a is weak at low temperature AF state, and should be sensitive to the 
occupation number of xy orbital. The less Uxy is, the stronger the transition is. Actually, this 
transition is much enhanced in the FM configuration as shown in Fig. 1(d). Therefore, with 
the increasing temperature, the enhancement of peak a can be understood by the increasing 
apical Ru-0 bond length and the increasing FM components in the spin configuration, both 
of which will suppress the xy occupation, while the reduction of peak (3 is mostly due to the 
spin disordering. The reduction of the gap at elevated temperature can be understood in 
terms of reduction in the AF long-range order. 
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B. Ca2_xSrxRu04 



Above results present firm evidence for tlie xy FO ordering for x=0.0. Now going from 
x=0.0 to 0.1, we see quantitative change of the electronic states (Fig. 1(e) and Fig. 2), while 
the ordering pattern does not change quahtatively. The elongation of the Ru-0 bond along 
the c-axis- (going from x=0.0 to 0.1) tends to suppress the occupation of xy orbital, and 
to destroy the AF order. As a result, we see from Fig. 2 the followings: 1) the energy gain 
of AF state with respect to the FM state reduces, being consistent with the reduced Neel 
temperature T/v observed experimentally; 2) the reduction of ordered moments; 3) electron 
transfer in the minority spin channel from the xy orbital (n^.j^) to the yzj zx orbitals (n^^/^^.)) 
4) redistribution of magnetization from the yzj zx orbitals [My^jz^ to the xy orbital iM^y). 
The reduction of n^^ is also visible from the calculated XAS spectra (shown in Fig.4(b)), 
where the relative intensity of two peaks (A and B) is just between those of Figs. 4(a) and 
(c). 

However for a;=0.5, where RuOe is quite elongated, and has only strong rotation around 
the c-axis without tilting, the situation is completely different (Fig. 1(f) and Fig. 2). These 
structural modifications produce basically two important changes in the electronic structure 
compared with the case of Ca2Ru04: 1) reduction in n\.y due to the elongation; 2) broadening 
of the yz/zx bands due to the absence of tilting. Nevertheless, the existence of RuOg 
rotation keeps the xy band narrow^^. This will contribute to the high density of states at 
the Fermi level of non-magnetic solution, leading to the Stoner-type FM instability. Since 
this instability mostly comes from the xy state, we observe from the FM solution (Fig. 1(f)) 
that the xy states are strongly spin polarized, opening a pseudo-gap, and contributing to the 
magnetization dominantly (Fig. 2). On the other hand, yz/zx bands are quite broad, and 
are located around the Fermi level with much reduced spin polarization compared with the 
case of x=0.0. In contrast to the 5* = 1 picture with the magnetic moment supported by the 
yz/zx orbitals for the x=Q case (Ca2Ru04), the present case may correspond to the 5'=l/2 
picture where the xy orbital contributes to the magnetic moment. This picture accounts 
well for the observations by the Curie- Weiss fitting with 5*= 1/2 of susceptibilit}*^ and the 
polarized neutron^ showing the spatial magnetic moment distribution of the xy character. 
The large spin polarization for the in-plane oxygen atoms and the negligibly small one for 
the apical oxygen atoms (Fig.2(c)) are the results coming from the xy orbital origin of the 
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Ru magnetic moment and are consistent with the experimenlj^. 

In contrast to the present study, the proposal in Ref.»lfl, which claims the contribution of 
yz/zx states to the magnetization, was obtained by neglecting the RuOg rotation, which is 
an important ingredient for obtaining our orbital-dependent picture. In their treatmenlfiS, an 
elaborate LDA+DMFT (dynamical mean field theory) scheme was used, but high symmetry 
structure of Sr2Ru04 was assumed. Then the effects of doping are simulated by increasing 
the size of U , instead of changing the band-width. However, we pointed out in our calculation 
that the strong RuOg rotation will reduce the band width of xy significantly, but not that of 
yz^ zx bands. This difference produces the main source of the discrepancy discussed above. 
The characteristic results of our picture for x=0.5 are the much suppressed (enhanced) peak 
A (B) in the XPS spectra (Fig.4(d)), and the structureless optical spectra below the charge- 
transfer peak (Fig.5(b)). These predictions should be confirmed by further experiments. 

Our picture for x=0.5 was obtained by the LSDA calculations which tends to overestimate 
the FM stability. Experimentally no FM long range order has been observed, while the 
significant enhancement of susceptibility down to 2K— clearly suggested the existence of 
strong FM correlation (at least for short range). Actually, the recent experiments down to 
O.SK"^ suggested the existence of weak FM component (cluster glass) at this doping. This 
may suggest that some magnetic orders may be energetically in near degereracy with FM 
state. At present we do not have any clear idea about what magnetic orders they may be 
and the problem is left for future studies. 

Finally, we discuss only briefly the case for the critical point x=0.2. lOK is already above 
the metal-insulator transition temperature^ in this case, and the system crystallizes in the 
L — Pbca phase with a long apical Ru-0 bond, suggesting the adjacency to the metallic side 
(x=0.5). However, in contrast to the case of a;=0.5, the tilting of RuOe still exists to reduce 
the width of yz/zx bands as well at x=0.2. This aspect tends to stabilize the insulating 
AF state. If we adopt LSDA, the former aspect is emphasized to make the case of x=0.2 
very similar to the case of x=0.5 as can be seen in Fig. 2. On the other hand, the LDA+U 
method with Ues =2.5 eV emphasizes the latter aspect to bring the AF state very close to 
the FM state with the energy difference given by -Efm ~ -f'AF = —3 meV. Note, however, 
that the similar LDA+t/ calculation for x=0.5 enhances the stability of the FM state with 
respect to the AF state and doubles the energy difference in Fig. 2(d). These calculations, 
though not conclusive, clearly suggest that the system with x=0.2 may be in the critical 
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situation. 



IV. CONCLUDING REMARKS 

We have presented a systematic picture for the orbital states in Ca2-xSr^Ru04 {0 < x < 
0.5). For the AF side {x <0.2), we conchide that the orbital ordering is of the FO type 
with the dominant xy occupation and that the yz/zx orbitals contribute to the magnetic 
moment. The electronic structure corresponds basically to 5"=! in the localized spin pic- 
ture. However, strong intersite hybridization due to extended Ru Ad orbitals significantly 
modifies the magnetic moment distribution and the present calculation agrees well with the 
experiments. For a;=0.5, our LSDA calculation predicts the system to be ferromagnetic and 
the narrow xy band with half filling may correspond to an S=l/2 localized spin picture. 
Experimentally the system shows strong tendency toward ferromagnetism but remains para- 
magnetic for a;=0.5 even down to very low temperature. Only below 0.3K, the existence of 
cluster glass phase with weak FM component was observed. The LSDA calculation seems 
to overestimate the stability of ferromagnetism in other cases also. This problem is left for 
future studies. However, the present result for the character of the spin-density distribution 
is consistent with the neutron data. We emphasized that the structural modifications such 
as rotation, tilting and fiattening of RuOg octahedron have to be fully taken into account 
in order to describe properly the systematic variation of the orbital states. Our results can 
explain other existing experiments, for example, XAS and optical conductivity, quite con- 
sistently. Furthermore, we made predictions for these experiments for the x—0.5 compound, 
which await experimental tests. 
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FIG. 1: The calculated PDOS of Ru-4d-t2g orbitals for various states. The solid red arrows 
indicate the main optical transition paths, while the dashed red arrow shows the transition which 
is significantly suppressed. Note the transitions indicated here mean that from the states at one 
Ru^^ site to the corresponding states at nearest-neighboring Ru'^"'" site. 
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FIG. 2: The calculated (a) orbital occupation numbers and (b) magnetic moments for the ground 
state of Ca2-3;Sr3;Ru04 with different doping x. M^y (My^/^^) denotes the magnetic moment 
associated with the xy (yz/zx) orbital and M^u-Mal = + 2.QMyzi^r^. The absolute value of 
magnetic moments from oxygens are shown in (c) (see text for further explanations). The panel 
(d) shows the total energy difference between the FM and AF states. The AF side (a;=0.0, 0.1) is 
connected with the FM side (x=0.2, 0.5) by dashed lines. 
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FIG. 3: The projected density of states of tetragonal Ca2Ru04 without the JT distortion in the 
NM state. The inset shows the calculated energy level splitting between the xy and the yz, zx 
states for CaRuOs and Ca2Ru04 with tetragonal Jahn- Teller distortion, which is defined as the 
bond length ratios (Ru-0 bond length along the c-axis versus that in the a6-plane). 
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FIG. 4: The calculated relative intensity of two XAS peaks (A and B located at 528.5 and 529.5 
eV experimentally) as a function of light incidence angle 9. 
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FIG. 5: The (a) experimental^ and (b) calculated optical conductivities. The interband transitions 
are calculated by using the Kubo formula. The simple average of FM and AF solutions is used to 
simulate the PM state at 295K. 
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